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基于 CO I 基因 序列 的 蜜蜂 属 系统 发 育 分 析 及 其 
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摘要 :【 目 的] 分 析 蜜 蜂 属 Apis 的 系统 发 育 关 系 , 并 在 此 基础 上 探讨 蜜蜂 属 钴 踊 方 向 、 姓 踏 声音 、 营 
集 环 境 、 梨 脾 结 构 的 祖先 状态 和 演变 过 程 。【 方 法 】PCR 扩 增 并 测定 中 国 分 布 的 蜜蜂 属 内 东方 蜜蜂 
A. cerana, 9 Z; Z% A. mellifera KZŽ A. dorsata, Z KZ% A. laboriosa, h Z% A. florea 和 黑 小 
蜜蜂 A. andreniformis 的 线粒体 CO 工 基 因 序 列 ,并 从 NCBI 数据 库 中 下 载 分 布 于 其 他 国家 或 地 区 的 
上 述 6 种 蜜蜂 以 及 绿 努 蜂 4，nulunsis、 苏 拉 威 西 蜂 4. nigrocinta 、 沙 巴 蜂 A. koschevnikovi, A X, X € 
XE A. dorsata binghami 的 CO 工 同 源 序 列 。 分 别 利用 最 大 简约 法 (MP)、 最 大 似 然 法 (ML) 和 贝 叶 斯 
分 析 ( BT) 依据 这 些 序 列 数据 构建 蜜蜂 属 系统 发 育 关 系 。 对 蜜蜂 属 上 述 种 的 三 蹈 语言 和 筑 业 行为 
进行 编码 并 作为 性 状 特征 标记 到 系统 发 育 树 中 ,利用 简约 法 对 祖先 状态 进行 追溯 。【 结果 】 系 统 发 
育 分 析 结 果 表 明 ,蜜蜂 属 可 划分 为 3 大 类 群 , 即 穴居 蜜蜂 类 群 (东方 蜜蜂 .西方 蜜 蜂 、 苏 拉 威 西峰 、 
绿 努 蜂 、 沙 巴 蜂 ) 、 大 蜜蜂 类 群 (大 蜜蜂 、 黑 大 蜜蜂 ) 和 小 蜜蜂 类 群 (小 蜜蜂 、 黑 小 蜜蜂 ) ;小 蜜蜂 类 群 
更 加 接近 于 祖 型 ,大 蜜蜂 类 群 和 穴居 蜜蜂 类 群 是 两 支 单 系 ; 炳 式 大 蜜蜂 是 独立 于 大 蜜蜂 和 黑 大 蜜蜂 
的 蜂 种 , 且 与 黑 大 蜜蜂 的 亲缘 关系 更 近 。 祖 先 状态 重建 结果 显示 :蜜蜂 属 祖 先 在 露天 环境 下 筑 造 重 
直 的 单 脾 , 且 在 传播 食物 或 梨 址 信息 时 束 水 平方 向 的 无 声 摆 尾 苗 , 有 喻 喻 声 的 舞蹈 及 筑 造 复 脾 是 后 
来 形成 的 。 【结论 ]CO 工 基因 可 作为 分 子 标记 用 于 分 析 蜜 峰 属 的 狂 蹈 和 筑 梨 行为 的 祖先 状态 及 进 
化 过 程 ;蜜蜂 筑 造 复 脾 、 跳 有 喻 喻 声 的 舞蹈 是 后 期 的 适应 性 进化 行为 。 

关键 词 : 蜜蜂 属 ; CO 工 基因 ; 系统 发 育 ; 舞蹈 语言 ; SLE] 
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Abstract: [ Aim] To analyze the phylogeny of the genus Apis, and to investigate the ancestral state and 


the evolutionary process of dance direction, dance sound, nesting environment and comb structure in Apis 
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( Hymenoptera ) 蜜蜂 总 科 ( Apoidea ) 蜜蜂 科 ( Apidae ) 
HREJ (Apini), 
昆虫 ,此 外 还 是 研究 社会 性 行为 的 关键 模式 生物 , 具 


based on the phylogenetic relationships. [Methods] The COI gene fragments of A. cerana, A. mellifera, 
A. dorsata, A. laboriosa, A. florea and A. andreniformis distributed in China were amplified by PCR and 
sequenced. The homologous sequences of COI gene of the above six species in other states or regions and 


A. 


The phylogenetic relationships were established using maximum parsimony ( MP), 


those of A. nulunsis , 
NCBI databases. 


maximum likelihood ( ML) and Bayesian inference ( BI) based on these sequence data. The dance 


nigrocinta, A. koschevnikovi and A. dorsata binghami were downloaded from 


language and nesting behavior of each Apis species mentioned above were coded and treated as trait 
markers, and added correspondingly into the phylogenetic trees. The ancestral state was traced with the 
ML model. [Results] The results of phylogenetic analysis demonstrated that Apis species can be divided 
into three clusters, i. e., cavity-nesting honeybees (A. cerana, A. mellifera, A. nigrocinta, A. nulunsis , 


and A. koschevnikovi ) , and A. 


florea, and A. andreniformis). The dwarf honeybees are more ancestral, and the giant and cavity-nesting 


giant honeybees (A. dorsata, laboriosa) and dwarf honeybees ( A. 


species are monophyletic. Our results also suggest that the giant bee 4. dorsata binghami is a separate 
species from A. dorsata and A. laboriosa and closer to A. laboriosa. Prototype inference indicated that 
ancestral honeybee species nested in the open with a vertical single comb. They disseminate information 
about food sources and nesting sites using a silent horizontal waggle dance. The behavior of dancing with 
buzz and nesting with multi-combs appeared subsequently. [ Conclusion] As a molecular marker, COI 
gene can be used to analyze the ancestral state and the evolution of dancing and nesting of Apis species. 
The buzzing dance language and multi-comb nesting behavior could be an outcome of adaptive evolution 
during the later stage. 


Key words: Apis; COI gene; phylogeny; dance language; nesting behavior 
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造 圆 形 或 椭圆 形 的 蜂 梨 , 梨 脾 项 部 用 来 贮存 花蜜 , 
较 平 坦 的 平台 (Barron and Plath, 2017; Beekman 
喜欢 在 中 空 的 
树干 或 洞穴 中 筑 造 竖 直 的 复 脾 (多 张 集 脾 ) ,但 在 某 
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有 重要 的 经 济 效益 .生态 效益 及 科研 价值 ( 曾 志 将 些 情况 下 ,西方 塞 蜂 也 会 依附 树枝 贷 造 露天 的 蜂 某 ， 
2007; Arenas et al., 2013; Silva et al., 2015) 目前 由 此 可 见 ,根据 环境 和 可 利用 的 机 会 不 同 ,蜜蜂 的 筑 
蜜蜂 属 内 种 类 分 为 3 大 类 群 , 即 穴居 蜜蜂 类 群 ( 东 — 梨 行 为 表现 出 较 大 的 灵活 性 (Engel and Schultz, 
方 蜜蜂 Apis cerana .西方 蜜蜂 Apis mellifera , 苏 拉 威 1997, Bhagavan et al., 2016) 。 近年 来 ,有 研究 者 相 
西 蜂 Apis nigrocinta , 绿 努 蜂 Apis nulunsis 沙巴 qe 继 根 据 蜜蜂 生物 学 信息 及 分 子 生 物 学 技术 ,对 祖先 
Apis koschevnikovi) .小 蜜蜂 类 群 ( 小 蜜蜂 Apis florea., 蜜蜂 是 在 腔 洞 中 筑 策 还 是 露天 营 策 进行 了 探讨 ,但 
黑 小 蜜蜂 Apis andreniformis) 和 大 蜜蜂 类 群 ( 大 蜜蜂 ”此 问题 仍然 存在 很 多 不 确定 性 (Raffiudin and 
Apis dorsata 、 黑 大 蜜蜂 Apis laboriosa.) ,而 大 蜜蜂 类 群 Crozier, 2007 ) 。 


中 的 炳 式 大 蜜蜂 Apis dorsata binghami 的 分 类 地 位 
目前 仍 存在 争议 (Alexander， 











在 社会 性 昆虫 中 ,蜜蜂 的 食物 信息 传播 非常 
杂 : 工 蜂 在 外 采集 到 较 好 的 蜜源 返回 蜂 集 后 ,会 








1991; Engel and 
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Schultz, 1997; Arias and Sheppard, 2005; Raffiudin 脾 上 做 出 有 规律 的 运动 ,以 鼓励 其 他 蜜蜂 出 Kee 
and Crozier, 2007) 。 蜜 蜂 属 内 种 类 具有 3 个 共同 的 ”或 告知 同伴 食物 的 位 置信 息 , 人 们 将 这 种 规律 性 运 
生物 学 特性 : 营 社 会 性 生活 ` 泌 蜡 筑 造 具有 六 角形 策 EKA PERI HC (Dyer, 2003; Toufailia, 2013) 。 
房 的 巢 脾 以 及 贮 密 积 极 ( 兽 志 将 ,2007 ) 。 蜜蜂 舞蹈 通常 发 生 在 垂直 于 地 面 和 水 平方 向 上 , 例 

大 蜜蜂 和 小 蜜蜂 类 和 群 的 蜂 梨 通常 是 露天 筑 造 的 ”如 :东方 蜜蜂 西方 蜜蜂 ,大 蜜蜂 和 黑 大 蜜蜂 的 梨 脾 






































PAR ,但 前 者 的 蜂 梨 通常 附着 在 树枝 或 崖 壁 下 LEER 
竖 直 垂下 ,顶部 没有 自由 空间 (Tan, 2007; Thakar 
and Tonapi, 2015) ;而 小 蜜蜂 类 群 一 般 围 绕 树 枝 筑 





















































是 竖 直 悬 挂 的 ,因此 它们 的 舞蹈 是 垂直 方向 的 ;而 小 
蜜蜂 和 黑 小 蜜蜂 的 舞蹈 发 生 于 梨 脾 上 方 的 平台 上 ， 
是 水 平方 向 的 (Beekman et al., 2008 ; Schürch et al., 
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2013) 。 有 研究 表明 ,将 西方 蜜蜂 的 梨 脾 水 平 放置 ， 
并 在 梨 脾 上 方 提供 一 个 光源 ,给 予 出 梨 的 工蜂 饲 喂 
奖励 ,工蜂 回 梨 后 跳 起 了 水 平方 向 的 舞蹈 ,因此 有 人 
认为 蜜蜂 舞蹈 的 方向 与 巢 脾 的 方向 紧密 相关 
( Benjamini et al., 2001 ) 。 人 们 针对 蜜蜂 “舞蹈 语 
言 ”进行 过 很 多 研究 ,但 通常 局 限于 破译 舞蹈 与 食 
物 源 之 间 的 关系 (Schiirch et al., 2013; Abrol, 2013; 
Price and Grüter, 2015) ,对 于 水 平 或 垂直 方向 的 舞 
蹈 究竟 哪 一 种 形式 是 祖先 的 特征 还 未 有 定论。 

蜜蜂 舞蹈 和 筑 梨 行为 的 演化 过 程 探究 需要 基于 
蜜蜂 属 的 系统 发 育 。 自 20 世纪 90 年 代 以 来 研究 者 
开始 从 形态 学 和 分 子 生物 学 方面 探讨 蜜蜂 属 种 类 间 
的 系统 进化 (Alexander，1991; Willis et al., 1992; 
Cameron et al., 1992; Engel and Schultz, 1997; 
Engel, 2000; Tanaka et al., 2001; 
Sheppard, 2005; 曹 联 飞 等 , 2012; Takahashi et al., 
2018) ,经 常 使 用 的 分 子 标记 有 线粒体 CO D, CO TI 
和 ND2 基因 ,但 Willis 4& ( 1992) 使 用 CO 本 基因 探 
究 蜜 蜂 属 部 分 物种 的 系统 发 育 关系 时 ,结果 显示 不 
文 持 现 有 沙巴 蜂 的 分 类 地 位 ,因此 他 认为 以 COI 基 
因为 分 子 标记 分 析 蜜 蜂 属 的 系统 发 育 具 有 一 定 的 局 
限 性 。CO I 基因 具有 长 度 适宜 、 进 化 速率 慢 、 富 合 
生物 系统 发 育 信号 等 优点 ,被 广泛 用 于 动 .植物 系统 
进化 研究 (Hebert et al., 2003) , Tanaka 等 (2001 ) 运 
用 CO 工 基因 研究 了 蜜蜂 属 内 5 种 蜜蜂 的 系统 发 育 
关系 ,研究 结果 支持 现 有 蜜蜂 的 分 类 地 位 。 

基于 上 述 情 况 , 本 研究 以 蜜蜂 属 9 种 蜜蜂 及 分 
类 地 位 存在 争议 的 炳 式 大 蜜蜂 为 研究 对 象 ,利用 
CO I 基因 重建 蜜蜂 属 的 系统 发 育 关 系 ,并 对 祖先 蜜 
蜂 的 舞蹈 形式 和 筑 梨 行为 进行 推断 ,分 析 蜜 蜂 舞 蹈 
语言 和 筑 梨 行为 的 演化 过 程 ,以 期 为 今后 蜜蜂 行为 
寺 征 的 进化 研究 提供 参考 ,并 为 蜜蜂 种 类 保护 及 野 
生 蜂 种 驯化 提供 依据 。 
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1 材料 与 方法 


1.1 实验 材料 与 DNA 数据 获取 

本 实验 于 2016 年 5 月 -2018 年 1 月 在 云南 农 
业 大 学 东方 蜜蜂 研究 所 分 子 实验 室内 完成 。 所 使 用 
的 蜜蜂 样本 分 别 来 源 于 云南 农业 大 学 东方 蜜蜂 研究 
所 实验 蜂 场 、 云 南 农 业 大 学 东方 蜜蜂 研究 所 标本 室 
和 云南 省 普洱 市 野外 ,包括 东方 蜜蜂 4，cerana、 西 
TRI A. mellifera 大 蜜蜂 4，dorsata 、 黑 大 蜜蜂 A. 
laboriosa、 小 蜜蜂 A. florea 和 小 蜜蜂 A. 















































andreniformis > 野外 采集 到 的 样本 全 部 浸泡 于 无 水 
乙醇 中 , 带 回 实验 室 后 置 于 -80Y 超低温 冰箱 中 保 
存 备用 。 

考虑 到 样本 采集 地 的 局 限 性 ,为 增加 每 种 蜜蜂 
的 生态 类 型 ,使 得 重建 的 系统 发 育 关系 更 具有 代表 
性 和 说 服 力 ,在 本 实验 采集 到 的 上 述 6 种 蜜蜂 样本 
中 每 种 每 个 采样 点 选择 1 头 用 于 测定 CO I 基因 ;再 
从 GenBank 数据 库 下 载 相关 的 蜜蜂 CO I 基因 序列 
用 于 比较 分 析 , 并 以 蜜蜂 科 麦 蜂 族 ( Meliponini) 麦 蜂 
Melipona scutellaris CO I 基因 (GenBank 登录 号 : 
KP202303) 作为 系统 发 育 分 析 的 外 群 。 本 实验 自 测 
和 下 载 序列 的 详细 信息 见 表 1。 
1.2 COI 基因 扩 增 及 序列 测定 

取 1.1 节 蜜 蜂 样 本 的 胸部 肌肉 ,利用 酚 / 氧 仿 法 
提取 基因 组 DNA ,用 TE 溶液 保存 提取 成 功 的 DNA, 
放置 于 -80% 超 低温 冰箱 中 保存 备用 (Zhao et al., 
2017), 
通过 PCR 扩 增 CO I 基因 片段 ,本 实验 中 使 用 
的 线粒体 CO I 基因 特异 性 扩 增 引物 为 LCOI490 
(forward ) : 5'-GGTCAACAAATCATAAAGATATTGG- 
3'; HC02198( reverse) : 5'-TAAACTTCAGGGTGACC 
AAAAAATCA-3' ( Tanaka et al., 2001) 。 由 生 工 生物 
工程 (上 海 ) 有 限 公 司 合成 后 ,10 000 r/min 离心 5 
min ,用 TE 溶液 溶解 至 终 浓 度 为 20 mmol/pL, 分 装 
后 放置 于 -80% 超低温 冰箱 中 保存 备用 。 

PCR 反应 体系 (40 uL): DNA 模板 200 ng, 10 
x 反应 缓冲 液 4 uL, dNTPs 4 uL, Taq 聚合 酶 2 U, 
上 下 游 引 物 (10 pmol/L) 各 2 uL, PCR 扩 增 程序 : 
95% 预 变 性 5 min; 94% 变性 20 s, 55% 退火 30 s, 
72% 延伸 30s, 30 个 循环 ; 72% 延伸 $ min, 终 止 反 
应 。PCR 产物 纯化 后 送 上 海 生 工 生 物 工程 有 限 公 
司 进行 双向 测序 。 
1.3 序列 比 对 及 系统 发 育 分 析 

测序 完成 后 ,首先 利用 SeqMan 程序 对 正 、 反 向 
CO 工 基因 序列 峰 图 进行 检查 校正、 修剪 并 生成 一 
致 序列 。 将 生成 的 一 致 序列 提交 到 NCBI 的 BLAST 
检索 ,与 已 提交 序列 进行 比 对 ,核实 是 否 为 目的 片 
段 ,之 后 将 所 测 CO I 基因 的 DNA 序列 上 传 至 
GenBank( 表 1) 。 联 合 GenBank 下 载 的 此 段 同 源 序 
列 ,将 1.2 T: PCR 扩 增 获得 的 以 及 从 GenBank 下 载 
的 26 条 CO 工 基 因 序 列 导入 MEGA6.0 中 ,在 默认 
参数 下 进行 Clustal W 计算 , 并 采用 最 大 简约 法 
( maximum parsimony , MP) 和 最 大 似 然 法 (maximum 
likelihood, ML) 构建 系统 发 育 树 , 系统 树 各 分 支 的 置 
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X1 COI 基因 序列 来 源 


Table 1 Sources of COI gene sequences 


种 名 样本 


Species 





Samples 


采集 地 
Collecting locality 


GenBank 登录 号 


GenBank accession no. 





东方 蜜蜂 Apis cerana Apis cerana 1 
Apis cerana 2 


Apis cerana 3 








方 蜜蜂 A. mellifera 








Apis mellifera 1 
Apis mellifera 2 
Apis mellifera 3 








大 蜜蜂 A. dorsata Apis dorsata 1 
Apis dorsata 2 
Apis dorsata 3 
炳 式 大 蜜蜂 A. dorsata binghami 


小 蜜蜂 A. florea 


Apis dorsata binghami 
Apis florea 1 
Apis florea 2 
Apis florea 3 


Phal 





BKCHEWE A. laboriosa Apis laboriosa 1 


Apis laboriosa 2 





黑 小 蜜蜂 A. andreniformis Apis andreniformis 1 
Apis andreniformis 2 
Apis andreniformis 3 


沙巴 蜂 A. koschevnikovi 





Apis koschevnikovi 1 
Apis koschevnikovi 2 
Apis koschevnikovi 3 
绿 努 峰 A. nulunsis Apis nulunsis 1 
Apis nulunsis 2 
Apis nulunsis 3 
苏 拉 威 西峰 A. nigrocinta 
麦 蜂 Melipona scutellaris 


Apis nigrocinta 


Melipona scutellaris 

















































































































中 国 云 南 Yunnan, China MH138082* 
泰国 Thailand AF153106” 
越南 Vietnam KM066955^ 

中 国 云南 Yunnan, China MH138083* 
俄罗斯 Russia MF100924^ 

摩洛哥 Morocco KY926883^ 

中 国 云南 Yunnan, China MH138081* 

马来西亚 Malaysia AY012725^ 
马来西亚 Malaysia DQ016080^ 
印度 尼 西 亚 Indonesia DQ020234^ 

中 国 云南 Yunnan, China MH138080* 
泰国 Thailand DQ016083^ 

中 国 云南 Yunnan, China KC170303^ 

中 国 云南 Yunnan, China MH138079* 
尼泊尔 Nepal AP018039” 

中 国 云南 Yunnan, China MH138078* 

未 知 Unknown EU163123^ 
印度 India KF736157"^ 
马来西亚 Malaysia AP017643” 
马来西亚 Malaysia AY754732? 
马来西亚 Malaysia DQ016101^ 
马来西亚 Malaysia MF565375^ 
马来西亚 Malaysia AF153099” 
马来西亚 Malaysia AP018157” 
印度 尼 西 亚 Indonesia KY799147^ 
未 知 Unknown KP202303^ 











a 本 实验 自 测 Obtained from this study; " H GenBank 下 载 Downloaded from GenBank. 


信 度 由 Bootstrap 1 000 次 重 抽 样 检 验 。 

使 用 MrBayes 3. 1. 2 进行 贝 叶 斯 系统 发 育 分 析 
(BD) ,在 分 析 前 使 用 MrModeltest2. 3 检测 序列 的 最 
适 碱 基 蔡 代 模 型 。BI 分 析 采 用 蒙特 卡 罗 马 尔 科 夫 
链 式 反应 (Markov Chain Monte Carlo, MCMC) 运行 
20 000 000 代 , 和 迭代 多 次 直至 BI 分 析 的 收敛 值 大 
0.01 ,每 1 000 代 取 样 1 次 ,将 最 初 的 10 000 个 样本 
作为 老化 样本 舍弃 ,而 将 剩余 的 10 000 个 样本 用 作 

















及 麦 蜂 M. scutellaris 的 行为 特征 的 祖先 状态 进行 重 


建 , 采 用 了 基于 分 支 树 追 溯 性 状 以 及 简约 法 祖先 状 
态 重建 。 


A 
TA 

















2: 结果 


2.1 蜜蜂 属 的 系统 发 育 
基于 本 实验 目测 和 下 载 自 GenBank 的 CO IÆ 





后 续 的 分 析 。 总 结 得 到 的 一 致 树 , 其 节点 的 可 靠 性 
采用 后 验 概率 (posterior probability，PP) 加 以 评估 
(Huelsenbeck et al., 2001 ) 。 
1.4 蜜蜂 行为 特征 祖先 状态 重建 

为 了 追溯 蜜蜂 属 舞 蹈 语言 和 筑 集 行为 的 进化 过 
程 ,首先 对 蜜蜂 属 舞 蹈 的 声音 、 方 向 及 营 集 环境 、 集 
脾 结 构 特 征 进行 赋值 ( 表 2)。 将 这 些 特征 作为 性 状 
标记 到 1.3 5 BI 分 析 获 得 的 系统 发 育 树 中 ,运用 
Mesquite3. 4. 0 程序 对 蜜蜂 属 9 种 蜜蜂 、 炳 式 大 蜜蜂 












































因 序 列 ( 表 1) ,运用 MEGA6.0( 使 用 最 大 简约 法 
(MP) 和 最 大 似 然 法 (ML) ) 及 MrBayes 3. 1.2 软件 ， 
以 麦 蜂 M. scutellaris 为 外 群 ,重建 了 蜜蜂 属 内 9 种 
蜜蜂 的 系统 进化 关系 ,结果 得 到 了 相似 的 拓扑 结构 ， 
都 支持 将 蜜蜂 属 划分 为 3 个 大 类 群 :穴居 蜜蜂 类 群 、 
大 蜜蜂 类 群 和 小 蜜蜂 类 群 ,其 中 小 蜜蜂 类 群 位 于 进 
化 树 的 基部 (图 1: A; 图 2) ， 说 明 其 最 接近 于 祖先 
类 型 ;大 蜜蜂 类 群 和 穴居 蜜蜂 类 群 为 后 来 分 化 形成 ， 
这 与 基于 形态 学 性 状 的 分 类 结果 一 致 (图 1 和 2)。 
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表 2 


蜜蜂 属 舞 蹈 语言 和 筑 梨 行为 特征 赋值 


Table 2 The coding of behavioral characters of dance language and nesting in the genus Apis 




























































































































































































特征 赋值 行为 特征 Behavioral characters 

Coding of EXE 75 TA] FAAN ERIE EXE 

characters Direction of dancing Sound during dance Nesting environment Nest structure 
1 不 跳舞 * No dancing 不 跳舞 * No dancing REK bie Open nesting 水 平 集 脾 和 * Horizontal combs 
2 水 平方 向 "" 无 声 的 舞蹈 洞穴 营 集 "… 位 于 支持 物 下 方 的 竖 直 单 脾 …i 

Horizontal direction Silent dance Cavity nesting One vertical comb beneath support 
3 垂直 于 地 面 ” 有 喻 喻 声 的 舞蹈" - 环绕 支持 物 的 竖 直 单 脾 … 
Vertical direction Dancing with buzz One vertical comb wrapped around support 
; ! - EFR Ft d i e 
Several vertical combs beneath support 

? 未 知 Unknown 未 知 Unknown 一 E 


: 无 性 状 No character. "Raffiudin and Crozier, 2007 ; "Engel and Schultz, 1997 ; 


, 2007; 


"Tan, 2017; 


*'Thakar and Tonapi, 2015; Ur vk 


l% , 2002; *Beekman et al., 2008 ; "Schürch et al., 2013; ' Barron and Plath, 2017; A 2013; Beekman and Oldroyd, 2018. 


A 


82 


76 一 一 大 蜜蜂 Apis dorsata 3 (DQ016080 





100 一 大 蜜蜂 Apis dorsata 1 (MH138081) 








大 蜜蜂 Apis dorsata 2 (AY012725) 




















o ”一 绿 努 峰 Apis nulunsis 3(AP018157) 
































97, —— 


E m 小 蜜蜂 Apis florea 1 (MH138080) 
小 蜜蜂 Apis florea 3 (KC170303) 





























82 


小 蜜蜂 Apis florea 2 (DQ016083) 








100 























84 





91 








94 


78 





74 小 蜜蜂 Apis florea 1 (MH138080) 
Æi Apis florea 2 (DQ016083) 


100 j 
小 蜜蜂 Apis florea 3(KC170303) 
RI 














LC 

















90 








100 








图 1 
Fig. 1 


bootstrap. 





























一 一 一 绿 努 蜂 Apis nulunsis 1 (MF565375) 
L— — — 东方 蜜蜂 Apis cerana 3 (KM066955) 
东方 蜜蜂 Apis cerana 1(MH138082) 

— 70 一 一 东方 蜜蜂 Apis cerana 2 (AF153106) 
苏 拉 威 西峰 Apis nigrocinta (KY 799147) 
方 蜜蜂 Apis mellifera 1 (MH138083) 
方 蜜蜂 Apis mellifera 2(MF100924) 
Jj Vi Apis mellifera 3(KY926883) 


黑 小 蜜蜂 Apis andreniformis 2 (EU163123) 
黑 小 蜜蜂 Apis andreniformis 1 (MH138078) 
88 黑 小 蜜蜂 Apis andreniformis 3 (KF736157) 
麦 峰 Melipona scutellaris (KP202303) 


黑 小 蜜蜂 Apis andreniformis 2 (EU163123) 
黑 小 蜜蜂 Apis andreniformis 1 (MH138078) 
85 黑 小 蜜蜂 Apis andreniformis 3 (KF736157) 
89 大 蜜蜂 Apis dorsata 3 (DQ016080) 
100 KÆ Apis dorsata | (MH138081) 
大 蜜蜂 Apis dorsata 2 (AY012725) 
炳 式 大 蜜蜂 Apis dorsata binghami (DQ020234) Giant honeybee 


"wl (—1 黑 大 蜜蜂 Apis laboriosa 1 (MH138079) 
100 黑 大 蜜蜂 Apis laboriosa 2 (AP018039) 


85 TRE Apis nulunsis 2 (AF153099) 
p OL SB Apis nulunsis 3 (AP018157) 
绿 努 蜂 Apis nulunsis 1 (MF565375) 


东方 蜜蜂 Apis cerana 3 (KM066955) 
东方 蜜蜂 Apis cerana 1 (MH138082) 
东方 蜜蜂 Apis cerana 2 (AF153106) 
沙巴 蜂 Apis koschevnikovi 2 (AY754732) 
沙巴 蜂 Apis koschevnikovi 1 (MF565375) 
沙巴 蜂 Apis koschevnikovi 3 DQ016101) 
苏 拉 威 西 蜂 Apis nigrocinta (KY799147) 
西方 蜜蜂 Apis mellifera | (MH138081) 
西方 蜜蜂 Apis mellifera 2 (MF100924) 
97 西方 蜜蜂 Apis mellifera 3 (KY926883) 
ZEE Melipona scutellaris (KP202303) 








大 蜜蜂 类 群 


HRK EH Apis dorsata binghami (DQ020234) Giant honeybee 
一 一 黑 大 蜜蜂 Apis laboriosa 1 (MH138079) 
100 黑 大 蜜蜂 Apis laboriosa 2 (AP018039) 
65 一 一 沙巴 蜂 Apis koschevnikovi 3 (DQ01610D 
一 一 沙巴 蜂 Apis koschevnikovi 1 CAPO17643) 
一 一 一 一 沙巴 蜂 Apis koschevnikovi 2(AY 754732) 
86 [一 一 绿 努 蜂 Apis nulunsis 2 (AF153099) 


穴居 蜜蜂 类 群 
Cavity-nesting 
honeybee 


小 蜜蜂 类 群 
Dwarf honeybee 





外 群 Outgroup 





小 蜜蜂 类 群 
IDwarf honeybee 


大 蜜蜂 类 和 群 





穴居 蜜蜂 类 群 
Cavity-nesting 
honeybee 





外 群 Outgroup 


最 大 简约 法 和 最 大 似 然 法 构建 的 基于 CO T 基因 序列 的 蜜蜂 属 系统 发 育 树 
Phylogenetic trees of the genus Apis based on the COI gene sequences using maximum parsimony and maximum likelihood 


A; MP fij Parsimony tree; B: ML 树 Likelihood tree. 图 中 分 支 处 数字 表示 自 展 支持 率 Values along branches in the figure represent the parsimony 
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从 系统 进化 分 析 的 结果 可 以 看 出 ,在 穴居 蜜蜂 
类 群 中 ,东方 蜜蜂 、 绿 努 蜂 、 沙 巴 蜂 和 苏 拉 威 西 蜂 紧 
密 地 聚集 在 一 起 :在 MP 和 ML 分 析 中 东方 蜜蜂 与 
绿 努 蜂 聚 为 一 支 , 自 展 文 持 率 为 78% (图 1: A, B), 
而 在 BI 分 析 中 ,东方 蜜蜂 与 苏 拉 威 西峰 眼 为 一 
支 ,BI 后 验 概率 为 0.76( 图 2), 却 都 显示 为 可 靠 的 


分 组 。 


















































在 大 蜜蜂 类 群 中 ,MP 和 ML 分 析 显 示 炳 式 大 蜜 
蜂 与 黑 大 蜜蜂 的 亲缘 关系 更 近 , MP 分 析 中 两 者 聚 
类 的 自 展 支 持 率 较 高 ,为 85% ,显示 聚 类 可 靠 (图 1: 
A) ;但 在 ML 分 析 中 此 节点 的 自 展 支持 率 仅 为 38% 
(图 1: B); mM BI 分 析 显 示 , 烦 式 大 蜜蜂 是 独立 于 黑 
大 蜜蜂 和 大 蜜蜂 之 外 的 蜂 种 ,但 后 验 概率 仅 为 
0. 53 ,支持 率 较 低 ( 图 2) 。 







































































1.00 
0.85 


0.53 


1.00. 


0.91 
1.00 


0.63 0.93 
1.00 
0.89 1.00 


0.98 


1.00 


0.98 


1.00 


0.86 
0.84 


1.00 














图 2 Wn 


Zl% Melipona scutellaris (KP202303) 
大 蜜蜂 Apis dorsata 3 (DQ016080) 
大 蜜蜂 Apis dorsata 1 (MH13808D 
大 蜜蜂 Apis dorsata 2 (AY012725) 


L MEE Outgroup 


大 蜜蜂 类 群 


炳 式 大 蜜蜂 Apis dorsata binghami (DQ020234 |Giant honeybee 


大 蜜蜂 Apis laboriosa 1 (MH138079) 
大 蜜蜂 Apis laboriosa 2 (AP018039) 
巴 蜂 Apis koschevnikovi | (AP017643) 
巴 蜂 Apis koschevnikovi 3 (DQ01610D 
沙巴 蜂 Apis koschevnikovi 2 AY 7154732) 
苏 拉 威 西峰 Apis nigrocinta KKY799147) 
东方 蜜蜂 Apis cerana 1 (MH138082) 
东方 蜜蜂 Apis cerana 2(AF153106) 
东方 蜜蜂 Apis cerana 3 (KM066955) 

ZEB Apis nulunsis 1 (MF565375) 

绿 努 蜂 Apis nulunsis 2 (AF153099) 

绿 努 蜂 Apis nulunsis 3 CAPO18157) 

西方 蜜蜂 Apis mellifera 1(MH138083) 
西方 蜜蜂 Apis mellifera 2 (MF100924) 
西方 蜜蜂 Apis mellifera 3 (KKY926883) 

小 蜜蜂 Apis florea 3 (KC170303) 

小 蜜蜂 Apis florea 1 (MH138080) 

小 蜜蜂 Apis florea 2 (DQ016083) 

四 小 蜜蜂 Apis andreniformis | (MH138078) 
小 蜜蜂 Apis andreniformis 2 (EU163123 
小 蜜蜂 Apis andreniformis 3 (KF736157) 


Xu du 


















































zm dm Hg 














穴居 蜜蜂 类 群 


Cavity-nesting honeybee 


小 蜜蜂 类 群 
Dwarf honeybee 


斯 法 构建 的 基于 COL 基因 序列 的 蜜蜂 属 系统 发 育 树 








Fig. 2 Phylogenetic tree of the genus Apis based on the COI gene sequences using Bayesian inference 
图 中 分 支 处 数字 表示 BI 树 的 后 验 概率 Values along branches in the figure represent the Bayesian posterior probabilities. 


2.2 ”蜜蜂 行为 祖先 状态 重建 

由 于 BI 分 析 法 得 到 的 系统 发 育 树 中 ,有 45% 
的 节点 得 到 了 最 大 后 验 概率 ( PP — 1. 00) ,85% 的 节 
点 得 到 了 0. 80 以 上 的 后 验 概率 ;而 MP 和 ML 分 析 
中 ,获得 最 大 自 展 值 的 节点 仅 占 33. 3% , 自 展 支持 
率 大 于 80% 的 节点 占 比 均 为 75%。 因 此 ,基于 BI 
分 析 得 到 的 系统 发 育 树 更 加 可 靠 。 所 以 ,本 实验 将 
BI 树 运 用 于 蜜蜂 属 的 舞蹈 语言 及 筑 梨 行为 的 祖先 
状态 重建 。 
2.2.1 蜜蜂 舞蹈 语言 :舞蹈 方向 的 状态 重建 结果 表 
明 ,蜜蜂 族 和 麦 蜂 族 的 共同 祖先 中 ,可 能 既 有 依靠 气 
味 痕迹 传递 食物 梨 址 信息 ,也 有 依靠 舞蹈 传播 的 。 
























































蜜蜂 属 祖先 的 舞 踊 最 可 能 是 水 平方 向 的 ,但 在 穴居 
蜜蜂 和 大 蜜蜂 类 群 的 共同 祖先 处 ,舞蹈 方向 发 生 了 
第 一 次 演变 ,转变 成 垂直 方向 ,之 后 一 直 处 于 稳定 状 
态 (图 3: A). 

蜜蜂 舞 中 声音 的 祖先 状态 推断 结果 如 图 3(B) 
所 示 。 从 图 中 可 以 看 出 ,祖先 蜜蜂 的 舞蹈 最 可 能 是 
无 声 的 , 带 有 喻 叭 声 的 蜜蜂 舞蹈 属于 派生 行为 ,此 状 
态 的 演化 同样 是 在 发 生 于 穴居 蜜蜂 与 大 蜜蜂 类 群 的 
共同 祖先 ,由 无 声 舞蹈 转变 为 再 有 喻 喻 声 的 舞蹈 ,而 
黑 大 蜜蜂 的 舞蹈 声音 仍然 遵循 祖先 状态 。 
2.2.2 BEARIN: 图 4(A) 为 蜜蜂 营 集 环境 的 
祖先 状态 重建 结果 ,图 中 显示 ,蜜蜂 族 与 麦 蜂 族 的 共 
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A B 
口 不 跳舞 No dancing (1) 苏 拉 威 西峰 Apis nigrocinta (KY799147) 国 不 跳舞 No dancing (1) 
国 水 平方 向 Q2) 东方 蜜蜂 Apis cerana 1 (MH138082) Bl XP Silent (2) 
Hori RU 东方 蜜蜂 Apis cerana 2 (AF153106) BAE 
orizontal direction : - ESKE Buzz (3) 
ging G) 东方 蜜蜂 Apis cerana 3 (KM066955) 




















沙巴 蜂 Apis koschevnikovi 3 (DQ01610D BI Unknown 


沙巴 蜂 4pis koschevnikovi | (AP017643) 
沙巴 蜂 Apis koschevnikovi 2 (AY754732) 
ZR Apis nulunsis 3 (AP018157) 
绿 努 蜂 Apis nulunsis 1 (MF565375) 
绿 努 蜂 Apis nulunsis 2 (AF153099) 
西方 蜜蜂 Apis mellifera 2 (MF100924) 
西方 蜜蜂 Apis mellifera 3 (KY926883) 
西方 蜜蜂 Apis mellifera | (MH138083) 
炳 式 大 蜜蜂 Apis dorsata binghami (DQ020234) 
黑 大 蜜蜂 Apis laboriosa 1 (MH138079) 
黑 大 蜜蜂 Apis laboriosa 2 (AP018039) 
大 蜜蜂 Apis dorsata 2 (AY012725) 
大 蜜蜂 Apis dorsata 1 (MH138081) 
大 蜜蜂 Apis dorsata 3 (DQ016080) 
小 蜜蜂 Apis florea 3 (KC170303) 
小 蜜蜂 Apis florea 1 (MH138080) 
小 蜜蜂 Apis florea 2 (DQ016083) 
黑 小 蜜蜂 Apis andreniformis | (MH138078) 
黑 小 蜜蜂 Apis andreniformis 2 (EU163123) 
小 蜜蜂 Apis andreniformis 3 (KF736157) 
Zł% Melipona scutellaris (KP202303) 


图 3 基于 BI 分析 构建 的 蜜蜂 舞蹈 的 进化 图 
Fig. 3 Evolution of honeybees dance based on the Bayesian topology 


A: 舞蹈 方向 Direction of dancing; B: £EJÉHU f Sound during dance. 图 例 中 数字 代表 特征 赋值 。Numbers in the figure legend represent the 
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character state values. 













A B 
BlEX ES Open nesting (1) 苏 拉 威 西峰 Apis nigrocinta (KY799147) BIKE JR Horizontal combs (1) 
MEA Cavity nesting 2) 东方 蜜蜂 Apis cerana 1(MH138082) -—— RETAN (2) 

东方 蜜蜂 Apis cerana 2 (AF153106) One verti comp beneath mpport 
东方 蜜蜂 Apis cerana 3(KM066955) 国 环绕 支持 物 的 紧 直 单 脾 (3) 
沙巴 蜂 Apis koschevnikovi 3 (DQ01610) One vertical comb wrapped around support 
沙巴 蜂 Apis koschevnikovi 1 (AP017643) 国王 位 于 支持 物 下 方 的 竖 直 复 脾 (4) 


沙巴 蜂 Apis koschevnikovi 2 (AY754732) 
绿 努 蜂 Apis nulunsis 3 (AP018157) 
绿 努 蜂 Apis nulunsis | (MF565375) 
绿 努 蜂 Apis nulunsis 2 (AF153099) 
西方 蜜蜂 Apis mellifera 2 (MF100924) 
西方 蜜蜂 Apis mellifera 3 (KY926883) 
西方 蜜蜂 Apis mellifera 1 (MH138081) 
炳 式 大 蜜蜂 Apis dorsata binghami (DQ020234) 
黑 大 蜜蜂 Apis laboriosa 1 (MH138079) 
黑 大 蜜蜂 Apis laboriosa 2(AP018039) 
大 蜜蜂 Apis dorsata 2 (AY012725) 
大 蜜蜂 Apis dorsata 1 (MH138081) 
大 蜜蜂 Apis dorsata 3 (DQ016080) 
小 蜜蜂 Apis florea 3 (KC170303) 
小 蜜蜂 Apis florea 1 (MH138080) 
小 蜜蜂 Apis florea 2 (DQ016083) 
黑 小 蜜蜂 Apis andreniformis 1 (MH138078) 
黑 小 蜜蜂 Apis andreniformis 2 (EU163123) 
黑 小 蜜蜂 Apis andreniformis 3 (KF736157) 
Æ$% Melipona scutellaris (KP202303) 














Several vertical combs beneath support 


图 4 基于 BI 分 析 构 建 的 蜜蜂 筑 梨 行 为 的 进化 图 


Fig. 4 Evolution of honeybees nesting based on the 
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A: 露天 筑 梨 或 洞 


represent the character state values. 





R 


I DA 





LEE gee 


EEA RREA 


WIER A R 


P> Ache (65 


UR TRAC BS 7I NIETEZSTR 

B , Tfi Ka AE 2I I P9] I IC BEER 

天 营 策 也 有 洞穴 筑 策 的 。 穴 居 蜜 蜂 类 群 的 祖先 是 洞 
穴 筑 梨 的 ,露天 营 梨 的 行为 消失 。 

蜂 集 构造 祖先 状态 推断 结果 显示 ,蜜蜂 属 祖先 

筑 造 垂直 的 单 脾 ,垂直 复 脾 属于 蜜蜂 进化 过 程 中 的 
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Bayesian topology 


Open nesting or cavity nesting; B: 集 脾 的 结构 Comb structure. 图 例 中 数字 代表 特征 赋值 。Numbers in the figure legend 





派生 状态 ,在 穴居 蜜蜂 的 祖先 处 发 生 此 种 演化 (图 
4: B), 4 村 
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结构 未 发 生 过 逆转 变 。 


本 实验 以 CO [基因 部 分 片段 为 分 子 标记 ,对 蜜 


1460 昆虫 学 报 Acta Entomologica. Sinica 6135 














蜂 属 的 系统 发 育 关 系 进行 了 重建 ,MP 和 ML 分 析 得 
到 的 系统 发 育 树 ,其 拓扑 结构 与 曹 联 飞 等 (2012 ) 利 
用 CO TL fü VD2 等 基因 得 到 的 蜜蜂 属 系统 发 育 关系 
相似 ,但 在 沙巴 蜂 的 分 类 地 位 上 略 有 不 同 。 本 实验 
中 ,沙巴 蜂 与 东方 蜜蜂 的 亲缘 关系 比 苏 拉 威 西 蜂 与 
东方 蜜蜂 的 更 近 , 推 测 沙巴 蜂 从 东方 蜜蜂 中 分 化 的 





















































KARIS HE JE fE, 这 可 能 是 一 种 适应 性 进化 。 
小 蜜蜂 类 群 和 黑 大 蜜蜂 的 舞蹈 在 白天 进行 ,因此 它 
们 依靠 视觉 来 接受 舞蹈 传递 的 信息 即 可 ;而 大 蜜蜂 
和 穴居 蜜蜂 虽 多 数 在 白天 舞蹈 ,但 有 时 也 会 发 生 在 
夜晚 (Abrol, 2013) , 且 穴 居 蜜 蜂 的 蜂 集 内 通常 光线 
非常 暗 , 所 以 这 两 者 除 依靠 视觉 外 可 能 还 要 依靠 听 



























































时 间 比 苏 拉 威 西 蜂 晚 ,而 曹 联 飞 等 (2012 ) 的 结果 显 
示 苏 拉 威 西 蜂 的 分 化 晚 于 沙巴 蜂 。 推 测 产 生 这 一 差 
异 的 主要 原因 在 于 使 用 的 分 子 标记 不 同 ,CO T 基因 
是 目前 动物 系统 进化 研究 应 用 最 频繁 的 分 子 标记 之 
一 ,具有 相对 保守 却 富 含 生物 系统 发 育 信号 等 优点 ， 
用 于 研究 种 间 的 系统 演化 关系 具有 很 大 的 优越 性 
(Hebert et al., 2003) ,而 Willis 在 使 用 CO 本 基因 对 
蜜蜂 属 部 分 种 类 进行 系统 发 育 重建 时 ,发 现 由 于 
CO 了 基因 进化 速率 较 快 等 原因 ,并 不 适用 于 蜜蜂 属 
的 系统 遗传 学 研究 ( Willis et al., 1992)。 因 此 我 们 
认为 本 研究 重建 的 蜜蜂 属 系统 发 育 关系 更 加 可 靠 。 

目前 炳 式 大 蜜蜂 的 分 类 地 位 仍然 不 太 明 确 , 虽 
然 有 学 者 把 炳 式 大 蜜蜂 归 为 大 蜜蜂 的 一 个 亚 种 
( Hepburn and Radloff, 2011) ,但 在 本 实验 的 MP 和 
ML 分 析 中 , 炳 式 大 蜜蜂 与 黑 大 蜜蜂 聚 为 一 文 , 且 在 
MP 分 析 中 得 到 了 85% 的 自 展 支持 率 。BI 分 析 显 示 
炳 式 大 蜜蜂 是 与 大 蜜蜂 、 黑 大 蜜蜂 并 列 的 物种 ,但 是 
后 验 概率 较 低 (0.53)。 这 与 Lo 等 (2010) 的 研究 结 
果 一 致 ,其 利用 BI 分 析 探 究 大 蜜蜂 类 群 中 的 系统 发 
育 关 系 时 发 现 炳 式 大 蜜蜂 是 独立 于 大 蜜蜂 和 黑 大 蜜 
蜂 之 外 的 物种 , 且 与 黑 大 蜜蜂 的 亲缘 关系 更 近 , 后 验 
概率 为 0.68。 但 由 于 GenBank 中 COT 基因 的 此 有 自 
同 源 序列 较 少 , 炳 式 大 蜜蜂 的 分 类 地 位 仍 需 进一步 
探讨 。 

蜜蜂 的 摆 尾 舞 被 认为 是 真正 具有 和 象征 意义 的 通 
讯 方 式 , 因 为 这 种 舞蹈 将 食物 源 的 时 间 和 空间 信息 
都 进行 了 编码 , 即 摆 尾 的 时 间 代 表 食 物 的 远近 ,蜜蜂 
头 部 的 朝向 与 重力 线 所 成 的 角度 代表 食物 位 置 与 太 
阳 所 成 的 角度 (Abrol, 2013) 。 本 研究 的 结果 显示 ， 
蜜蜂 属 祖先 的 摆 尾 舞 最 可 能 发 生 在 水 平方 向 上 ,最 
接近 祖 型 的 小 蜜蜂 类 群 在 传递 食物 和 和 集 址 信息 时 ， 
常 在 集 脾 上 方 的 平台 上 作 摆 腹 直 线 跑 动 ,而 运动 的 
朝 问 直 接 指示 了 目标 的 方位 ( 周 冰峰 , 2002) ,这 相 
比 于 穴居 蜜蜂 和 大 蜜蜂 类 群 依赖 重力 来 传递 方向 信 
息 无 疑 是 非常 节省 能 量 的 。 但 随 着 蜜蜂 的 进化 , 摆 
尾 舞 由 水 平方 向 进化 为 垂直 方向 ,产生 这 一 进化 的 
原因 还 不 明确 , 仍 需 要 深入 研究 。 属 内 祖先 蜜蜂 的 
舞蹈 是 无 声 的 ,但 大 蜜蜂 与 穴居 蜜蜂 的 共同 祖先 却 












































































































































































































































觉 接受 舞蹈 信息 。 

六 居 蜜 蜂 类 群 倾向 于 在 树木 的 腔 洞 或 洞 入 中 筑 
造 复 脾 ,这 是 一 种 衍生 的 行为 ,因为 祖先 蜜蜂 是 在 开 
放 的 环境 中 贷 造 单 脾 , 这 一 结果 与 Raffiudin 和 
Crozier(2007) 的 研究 结果 一 致 。 比 较 原 始 的 小 蜜蜂 
类 群 和 大 蜜蜂 类 群 分 布 于 亚洲 热带 地 区 ,分 布地 的 
年 平均 气温 在 12 以 上 , 整 群 蜜蜂 只 有 一 片 巢 脾 ， 
以 蜜蜂 黎 盖 集 脾 的 方式 就 能 维持 梨 温 ( 兽 志 将 ， 
2007)。 随 着 蜜蜂 分 布地 的 扩张 ,有 些 种 群 开始 登 
陆 温 带 地 区 ,此 地 区 冬季 的 温度 较 低 ,对 露天 营 梨 的 
蜜蜂 生存 提出 了 较 大 的 挑战 。 因 此 在 洞穴 内 筑 造 复 
脾 是 蜜蜂 种 群 的 一 大 进步 ,由 若干 片 梨 脾 组 成 的 复 
脾 不 仅 使 蜂 集 内 部 环境 更 加 稳定 , 且 避 人 免 阳光 直射 
对 蜂 儿 的 发 育 也 十 分 有 利 ;洞穴 筑 集 使 得 蜂 群 与 外 
界 的 接触 面积 减 小 ,在 很 大 程度 上 减少 了 外 界 环境 
(包括 天 敌 在 内 ) 的 干扰 ,并 且 随 着 防卫 蜂 的 数量 减 
D ,采集 蜂 的 数目 增加 ,为 蜜蜂 完整 社会 性 的 形成 也 
提供 了 可 能 ( 曾 志 将 , 2007) 。 
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